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In this letter we present, in a number conserving framework, a model of interacting fermions in a 
two-wire geometry supporting non-local zero-energy Majorana-like edge excitations. The model has 
an exactly solvable line, on varying the density of fermions, described by a topologically non-trivial 
ground state wave-function. Away from the exactly solvable line we study the system by means of the 
numerical density matrix renormalization group. We characterize its topological properties through 
the explicit calculation of a degenerate entanglement spectrum and of the braiding operators which 
are exponentially localized at the edges. Furthermore, we establish the presence of a gap in its single 
particle spectrum while the Hamiltonian is gapless, and compute the correlations between the edge 
modes as well as the superfluid correlations. The topological phase covers a sizeable portion of the 
phase diagram, the solvable line being one of its boundaries. 


Introduction — Large part of the enormous attention 
devoted in the last years to topological superconduct¬ 
ors owes to the exotic quasiparticles such as Majorana 
modes, which localize at their boundaries (edges, vor¬ 
tices, ... ) HIE] and play a key role in several robust 
quantum information protocols |3]. Kitaev’s p-wave su¬ 
perconducting quantum wire |3j provides a minimal set¬ 
ting showcasing all the key aspects of topological states 
of matter in fermionic systems. The existence of a so- 
called “sweet point” supporting an exact and easy-to- 
handle analytical solution puts this model at the heart 
of our understanding of systems supporting Majorana 
modes. Various implementations in solid state Eli] and 
ultracold atoms 1315 ] via proximity to superconducting 
or superfluid reservoirs have been proposed, and experi¬ 
mental signatures of edge modes were reported |9|. 

Kitaev’s model is an effective mean-field model and its 
Hamiltonian does not commute with the particle number 
operator. Considerable activity has been devoted to un¬ 
derstanding models supporting Majorana edge modes in 
a number-conserving setting [lOHIl] , as in various exper¬ 
imental platforms (e.g. solid state Hnmi] or ultracold 
atoms [HIS]) this property is naturally present. It was 
realised that a simple way to promote particle number 
conservation to a symmetry of the model, while keeping 
the edge state physics intact, was to consider at least 
two coupled wires rather than a single one [Mil]. How¬ 
ever, since attractive interactions are pivotal to gener¬ 
ate superconducting order in the canonical ensemble, one 
usually faces a complex interacting many-body problem. 
Therefore, approximations such as bosonization |inH12| . 
or numerical approaches m were invoked. An exactly 
solvable model of a topological superconductor in a num¬ 
ber conserving setting, which would directly complement 
Kitaev’s scenario, is missing (see however |14|). 

In this letter we present an exactly solvable model 
of a topological superconductor which supports exotic 


Majorana-like quasiparticles at its ends and retains the 
fermionic number as a well-defined quantum number. 
The construction of the Hamiltonian with local two-body 
interactions and of its ground state draws inspiration 
from ideas on dissipative state preparation for ultracold 
atomic fermions here applied to spinless fermi¬ 

ons in a two-wire geometry. The solution entails explicit 
ground state wave-functions, which feature all the main 
qualitative properties highlighted so far in approximate 
analytical Hanani] and numerical HSUS] studies for 
this class of models, with the advantage of being easy-to- 
handle. 

In particular, we establish the following key features: 
i) The existence of one/two degenerate ground states de¬ 
pending on the periodic/open boundaries with a two-fold 
degenerate entanglement spectrum; ii) the presence of 
exponentially localized, symmetry-protected edge states 
and braiding matrices associated to this degeneracy; iii) 
exponential decay of the fermionic single particle correl¬ 
ations, even if the Hamiltonian is gapless with collective, 
quadratically dispersing bosonic modes; iv) p-wave super¬ 
conducting correlations which saturate at large distance. 

By tuning the ratio of interaction vs. kinetic energy 
of our model, we can explore its properties outside the 
exactly-solvable line. The full phase diagram (Fig. 
is obtained by means of density matrix renormalization 
group (DMRG) calculations. The exactly solvable line is 
found to stand between a stable topological phase and a 
phase-separated state. 

The model — We begin by recapitulating some prop¬ 
erties of the Kitaev chain, whose Hamiltonian reads |1] 

Hk = ^ [ — JOjOj+i — Adjdj+i H.c. — — 1/2)]. 
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Here, J > 0 denotes the hopping amplitude, /i and A the 
chemical potential and the superconducting gap, respect¬ 
ively; are fermionic annihilation (creation) operators 
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on site j, and fij = ajoj. This model has i) two density- 
driven phase transitions from finite densities to the empty 
and full states at |/r| = 2 J (for A 7 ^ 0), and ii) a transition 
driven by the competition of kinetic and interaction en¬ 
ergy (responsible for pairing) at A/J = 0 (for |^| < 2J). 
For 1^1 < 2J and A 7 ^ 0, the ground state is unique for 
periodic boundary conditions, but twofold degenerate for 
open geometry, hosting localized zero-energy Majorana 
modes. This topological phase is symmetry protected by 
total fermionic parity P = (—1)^, where N = 

Let us focus on the so-called “sweet point”, namely /r = 
0, and A = J > 0 and real, which enjoys the property 
Hk = = “j + ^]+i 

and Aj = dj — dj+i {II is defined identifying L -\- 1 = 
1). For open geometry, the two ground states with L 
sites satisfy = 0 , for 1 < j < L — 1 , and can be 

written [SD] as the equal weighted superposition of all 
even (e) or odd (o) particle number states: 


/ e(o) 


i: 


J2n(2n+1) 


( 1 ) 


{j2n(2n + l)} 


Here |jm) = |vac) with ji < j,+i {ji = 

1, . . . ,1/) and Me,L — J2n ( 271 )’ Afo,L = J2n (2n+l) ’ 

We now turn to a number conserving version of this 
model on a single wire m- Indeed, the following model 
reduces precisely to the above scenario upon perform¬ 
ing a naive BCS mean field treatment. Consider the 
Hamiltonian with Lj = C^Aj, whose 

exact ground state wave-functions can be obtained as 
follows. Since Aj|V’)e(o) = -C']l^)e(o), l'0)e(o) are also 
ground states of Li\tlj)e(o) = 0 because (Cj)^ = 0 . 
As Li conserves the particle number, [Li,N] = 0, we 
can classify ground states for each fixed particle num¬ 
ber sector N by number projection, A) = Pw^eio)- 
This is implemented by choosing the state with 2n = N 
(or 2n -I- 1 = N) in the sum over n in Eq. Q, and ad¬ 
justing the normalization to Nl,n = (^)- The positive 
semi-definiteness of implies that these states, hav¬ 
ing zero energy eigenvalue, are ground states. However, 
once N is fixed, the ground state |'I', N) is unique, as fol¬ 
lows from the Jordan-Wigner mapping to the Heisenberg 
model [21]. The topological twofold degeneracy is lost. 

Guided by the previous analysis, we construct an ex¬ 
actly solvable topological two-wire model with fermionic 
operators d^P, Up. In addition to those involving each 
wire j we introduce new operators 

= CljAbj + CljAaj- The Hamiltonian 

E (2) 

cic—a,b,I j — 1 

coincides with the A = 1 point of the following more 
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Figure 1. Phase diagram for the number and local parity 
conserving two-wire model (§ as a function of A and filling 
V = N/2L obtained through DRMG simulations. The exactly 
solvable topological line is at A = 1 (another, trivially solvable 
non-topological line is at A = 0). For A > 1, the system 
undergoes phase separation (see the density profile (ri“) in 
the inset). For 0 < A < 1 and u ^ 0, 1, the system is 
in a homogeneous topological phase (see inset). The phase 
diagram is symmetric with respect to half filling = 1/2 due 
to particle-hole symmetry of H\. 


general model: 


L-l 

= (ajdj+i-fH.c.) -(h“ -f h“+i) -bAh“h“+;^ 

j — l,06—a,b 

L-l _ 

(n“ -b «“+i)(n5 + - ipaj+ib]bj+i 


- 2 AE[(^ 

i=i 

Ha. a 


+ d]dj+ibp^bj - 2b]b]^j^dj+idj + H.c.) . (3) 


H\ conserves the total particle number N = Na + Nb 
and the local wire parities Pa,b = (— 1 )'^“’*’, which act 
as protecting symmetries for the topological phase. The 
coupling A tunes the relative strength of the kinetic and 
interaction terms similarly to A/J in JJk- Although only 
A = 1 is exactly solvable, we will later consider A 7 ^ 1 
to explore the robustness of the analytical results. The 
phase diagram is anticipated in Fig. 

Exact results for A = 1 — For a fixed particle num¬ 
ber N and open boundaries, the ground state of H is 
twofold degenerate, due to the freedom in choosing the 
local parity. For even N, the ground states read 


N/2 

\4’l{N)) ee—■^eepN'^^ E \32n) a ® Wn- 2n) b : (4) 

{QN-2n} 

Ar/2-1 

\Ul{.N))^^= Nppm ^ ^ \j2n+l)a ® \qN-2n-l)b 

"=0 {j 2 „+l}. 

{^lV-2n-l} 






















































































































































3 


where Afee,L,N — Sn=0 ( 2 ™) (Ar-2n)’ ■^oo,L,N — 
T.n=o~^ { 2 n+i){N-L-i)- ^he states \j)a and are 
simple generalizations of the states |j) defined in Eq.Q 
to the wire a and b respectively. These describe the cases 
of even (ee) or odd ( 00 ) particle numbers in each of the 
wires. For odd N, the ground states |'0L(.^))eo(oe) ’'^ith 
an even (odd) number of particles in either wire take the 
identical sum structure as above with the normalization 
36. in both cases. The wave-functions Q are the 
unique ground states of the model ^3]. An interesting 
interpretation of \'4’l{N)) is in terms of number pro¬ 
jection of the ground state of two decoupled even-parity 
Kitaev chains |G) = |'0)g 0 |'*/’)e- 

HdN)).e ^ Pn |G) ; \i’Lmoo « |G); 

\MN)),,^PnIV\G); (5) 

where and are the zero-energy modes of the de¬ 
coupled Kitaev wires at half filling. This interpretation 
provides intuition that the two-fold ground-state degen¬ 
eracy is absent for periodic boundary conditions: since 
on a circle i?K has a unique ground state, the ground 
state of H with N particles is also unique |23| . 

Important evidence of a topologically nontrivial bulk 
state is obtained from the double degeneracy of the 
entanglement spectrum, which we now compute for 
one of the wave-functions Q. To this end, we con¬ 
sider the reduced state of I sites on each wire pi = 
Tr(L_i) [ {Pl{N)\^^]. Taking the symmetries 

into account, it can be written in diagonal form as |23| 

min{2l,N) 

Pl= {Pli.Nl)\aa' (6) 

TV;—0 <7^a' 

with the following nonzero eigenvalues: for Ni 

(oo),l,NiJ^ee(oo),L-l,N-Ni /Mee,L 

whereas for Ni odd 

the odd-particle number sector the entanglement spec¬ 
trum is manifestly twofold degenerate. In the even one, 
such degeneracy appears in the thermo^namic limit: 

X(e6l]v,)/X(ooIk) ^ 1 E3] and Fig. ^). 

An interesting insight is provideci by Oj = 
{tpL{N)\^^d^jbj |'0L(fV))ee, where is the only single¬ 
site operator which commutes with N and changes the 
local parities Pa,b, so that the two ground states can be 
locally distinguished. The calculation of such matrix ele¬ 
ments leads to a lengthy combinatorial expression [53] 
and is shown in Fig. H- We interpret the exponential 
decay of Oj into the bulk as a clear signature of local¬ 
ized edge modes with support in this region only. At half 
filling the edge states are maximally localized, but away 
from half filling the number projection increases the loc¬ 
alization length. In the thermodynamic limit, this length 
diverges for i/ = N/2L —)■ 0,1, indicating a topological 



Figure 2. Analysis of model H. (a) Entanglement spectrum 
for a reduced state pi with I = 60 for L = 240. (b) DMRG 
results for the scaling of the gap computed at fixed parity, 
which is compatible with L~^ (dashed line); here the num¬ 
ber of kept states is m = 400. (c) Localization of the edge 
mode computed via | (V’n(-X)loo IV'i(A')),,^ |. (d) Single¬ 
fermion edge correlations |(a]aj)| computed for a system of 
size L = 240. The wave-function is shift invariant, such that 
\{dldi+j)\ = Ka]aj)| {i + j< L). 


phase transition. We emphasize that this exponential be¬ 
havior is different from naiiii, reporting algebraic loc¬ 
alization of the edge states, but similar to [T31[IS|. Non¬ 
local correlations of edge states are another clear indic¬ 
ation of topological order and can be proven via {a\dj), 
which is sizeable both at j ^ 1 and j ^ L (see the ana¬ 
lytical expression in [53] and Fig. |^). 

Furthermore, the Hamiltonian is gapless and hosts long 
wavelength collective bosonic excitations, while the single 
fermion excitations experience a finite gap. This is a cru¬ 
cial property of the ground state; the absence of gapless 
fermion modes in the bulk ensures the robustness of the 
zero energy edge modes, in analogy to non-interacting 
topologically non-trivial systems. The gapped nature of 
single fermion excitations is established via the exponen¬ 
tial decay of the fermionic two-point function, e.g. {d\dj). 
Again, the resulting formula is a lengthy combinatorial 
expression [23] , evaluated numerically for very large sizes 
and plotted in Fig. Hi- For V 0,1, the correlation 
length diverges, indicating the vanishing of the fermion 
gap and a thermodynamic, density-driven phase trans¬ 
ition in full analogy to the Kitaev chain. 

On the other hand, the analysis of the superfiuid cor¬ 
relations demonstrates the existence of gapless modes. 
The p-wave nature of these correlations follows from the 
correlation of the pairing operator djj.idj. A direct cal¬ 
culation [53] shows a saturation at large distance 


\ L^oo 2/-I \2 

aK+iaj+iaj) —^ V [I- v) 


( 7 ) 
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Figure 3. DMRG results for model Hx- (a) Entanglement 
spectrum for a reduced state pi with I = 100 for L = 240 
(m = 300). (b) Algebraic scaling of the gap computed at 
fixed parity, which is compatible with L~^ (dashed line). Here 
m = 420. (c, d) Single-particle |(djaj)| and p-wave superfluid 
correlations at distance r = \i—j\ computed in 
the bulk of system with L = 240 (m = 300). Analogous data 
were obtained for other values of v, A (red circles in Fig. [^. 


Similar expressions hold for cross-correlations between 
the wires. The finite asymptotic value in Eq. Q hints 
at the absence of bosonic modes with linear dispersion, 
which would lead to algebraic decay. A DMRG analysis 
of the excitation spectrum of H for system sizes up to 
L = 144 demonstrates a vanishing of its gap ~ L~^ 
(Fig#). This indicates the presence of collective excit¬ 
ations with quadratic dispersion. Further support to this 
statement is provided from the fact that Eq. © without 
the wire coupling term reduces to the XXZ model at the 
border of its ferromagnetic phase, which hosts quadratic- 
ally dispersing spin waves, uj q^. This dispersion, with 
dynamic exponent z = 2, gives rise to an effective phase 
space dimension des = z + 1 = 3 zero temperature, ex¬ 
plaining the constancy of superfluid correlations due to 
the absence of a divergence in the soft mode correlators. 
This finding is special for A = 1. 

Non-abelian statistics — We now proceed to demon¬ 
strate that the edge modes obey a non-abelian stat¬ 
istics completely equivalent to that of Majorana fer¬ 
mions “ i.e., Ising anyons. Consider the operator 
BaR,bR{j) = + ZaR^bR,j)lV2 with j < L/2, where 

ZaR,bRj = (Z)p=l[n5=0 ZaR,bR,q]XaR,bR,p)/B{j), with 
XaR,bR,j = ~ yaR,bR,j = 

+ (1 - <+i_j)(l - <+i-j) for j > 
0, YaR,bR,0 = i and X, = - [^2# (1 _ 

BaR,bR{j) is thus exponentially localized at the right edge 
of the ladder and an analogous operator BaL,bL{j) can be 
defined at the left edge through the transformation map¬ 


ping an operator at site I to site L + l — l (and viceversa). 
Similarly, the operators BaR^ahij) and BbL,bR{j) can 
be defined through the transformations —)■ —iai 

and —>■ —ibi, respectively. In general, one can 

define operators BmA,m'A'ij) with m,m' = a,b and 
A, A' = L, R. These operators have the following key 
properties: They (i) are exponentially localized at the 
edges, (ii) act unitarily in the ground state subspace, 
(iii) are particle number conserving, and (iv) most im¬ 
portantly, provide a representation of Majorana braid¬ 
ing operators. From this we conclude that the local¬ 
ized edge modes behave as non-abelian Majorana fer¬ 
mions [32], in full analogy to the case of two neighbor¬ 
ing Kitaev wires. Properties (i)-(iii) are demonstrated 
in |33| , whereas here we focus on (iv). Strictly speaking, 
properties (ii) and (iv) are only true apart from an error 
which is exponentially small in j and L, which can al¬ 
ways be made negligible. In this case we can define the 
braiding operator = BmA,m'A'{j)- We initialize 

the system in the state |V’L(.^))ee then perform two 
braiding operations on the edges in different sequences. 
If we consider for example RaR,aL and RaR.bR we ob¬ 
tain [RaR,aL, BaRAR] |V'L(^))ee = ^ I'0 l(^))oo E3| which 
demonstrates the non-abelian character of these opera¬ 
tions. Moreover, this is the pattern that the conventional 
braiding operators produce on two neighboring Kitaev 
wires R'^A,m'A' = eTT'™A7„'A' = (I -b 7mA7m'A')/v^. 
where 7mA are Majorana operators fulfilling the Clif¬ 
ford algebra appearing at the left and right (A = L, R) 
edges of two Kitaev wires m = a,b. This pattern co¬ 
incides for the application of [i?m,A,m'A', RnX.n'X'] on all 
\'4’l{N ))states (see e.g. jH]). In other words, the op¬ 
erators RmA,m'A' form a number-conserving representa¬ 
tion of Majorana braiding operators on the ground state 
subspace. 

Numerical results — To further explore the status of 
these results, we now move to the full model Hx away 
from the solvable line A = 1. The study is performed 
with DMRG on systems with sizes up to L = 240 and 
open boundary conditions. 

We first establish the absence of a topological phase for 
A > 1. The density profile, shown in the inset of Fig.[^for 

= 0.5 and A = 1.01, displays a clear phase-separation 
tendency. Analogous data are obtained for other val¬ 
ues of iz (see dark crosses in Fig. [^. These results can 
be intuitively understood considering that i?A>i without 
interwire coupling can be mapped to a gapped ferromag¬ 
netic XXZ model with domain walls dual to fermionic 
phase separation. 

For A < 1, simulations support the existence of a homo¬ 
geneous phase (Fig.[^. Note that A = 0 is a free-fermion 
point trivially non-topological. For A 0 we observe: i) 
two quasi-degenerate ground states with different relat¬ 
ive parity and same particle numbers, ii) degenerate en¬ 
tanglement spectrum, iii) a gap closing as L~^ for fixed 
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parity, iv) exponentially decaying single-fermion correla¬ 
tions, v) power-law decaying superfluid correlators. Plots 
in Fig. 1^ display our numerical results. Simulations at 
lower Ailing ^ —>■ 0 and small A are more demanding, 
owing to the increasing correlation length of the system. 
The numerics is consistent with the phase diagram in 
Fig. [2 exhibiting a topological phase delimited by three 
trivial lines at A = 0, = 0 and ly = 1 and an inhomo¬ 

geneous non-topological phase for A > 1. The exactly 
solvable topological line at A = 1 serves as a boundary; 
the flxed-z/ phase diagram is reminiscent of the ferromag¬ 
netic transition in the XXZ model. 

Conclusions — We presented an exactly solvable 
two-wire fermionic model which conserves the number of 
particles and features Majorana-like exotic quasiparticles 
at the edges. Our results can be a valuable guideline 
to understand topological edge states in number con¬ 
serving systems. For example, the replacement hi —>■ 
results in a one-dimensional spinful Hub¬ 
bard Hamiltonian without continuous spin rotation, but 
time reversal symmetry. The resulting model with an 
exactly solvable line belongs to the class of time reversal 
invariant topological superconductors |25| . analyzed in a 
number conserving setting recently m, with edge modes 
protected by the latter symmetry. Moreover, exactly 
solvable number conserving models can be constructed 
in arbitrary dimension. 

During the Anal step of preparation, we became aware 
of similar results obtained by Lang and Biichler |26| . 
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Supplemental Material to 

Localized Majorana-like modes in a number conserving setting: 

An exactly solvable model 

Fernando lemini, Leonardo Mazza, Davide Rossini, Rosario Fazio, Sebastian Diehl 


In this Supplemental Material we provide additional information about some details of the analytical results for 
the exactly solvable two-wire topological system which have been omitted from the main text. 


TWO-WIRE GROUND STATE 


Hard-wall boundary conditions 


In this section we show that the wave-functions |'</’L(.^))ee(oo) @ main text are the only ground states 

of the two-wire Hamiltonian H\^i in the presence of hard-wall boundary conditions. Our proof actively constructs all 
of the zero-energy eigenstates of the Hamiltonian, which are the lowest-energy states because Hx-i > 0. Such states 
are obtained projecting the grand-canonical ground state of two decoupled Kitaev chains onto a given particle-number 
sector. 

Let us first consider only the operators {LojA&j} the corresponding parent Hamiltonian Hab = 

J2a=a b ’"^hich Corresponds to two decoupled wires. We know that the ground states of each wire are 

given by lb):- Hamiltonian Hab thus has a ground space spanned by 

{p:ph-nW>wiY ^ (SI) 

l. J n—O 

and (cr, cr') take the values (e, e) and (o, o) when N is even and (e, o) and (o, e) when N is odd. 

An important relation holds because |^/>)“ 0 = 0- Upon the insertion of the identity operator J2n n'=o PnPn' 

we get 


E 

n,n' 


(Cl .p: 




n—1 


-1-4 P 


)Pn' 


= 0 . 


(S2) 


Each of the elements in the above sum must vanish independently because of orthogonality, and the important relation 


/Aft pa p6 


, = lb)“ o , Vn, n' 


(S3) 


is derived (the same holds for the wire b). 

Let us now compute the most generic iV-fermions state such that = 0. In general, \(j)N) must be in the 

kernel of Hab'- 


N 


\4'n) — El El Xn,tT,iT'PnPN-r. 


F E 


^n,cr,(T' I 


= 1 . 


(S4) 


n—0 c 


Imposing now that Hi \(I)n) = = 0, we obtain that: 

N 

LL3 \M = (CljAb,, + ClAa,,) E E ^n,a,a'P;iP^N-n \P)l ® \P)l , 

n=0 <7,(7' 

= CPAb,, - CPAb,, E E ^n,a,a'P;i-2pN-n+2 Wl « 

n (7,(7' 

\b i./.\a ^ \_i.\b 


)E 


“ P'ljAbJ El E1(^".^.P’' Xn+2,a,a')PnPN-r 

n (7,(7' 


, =0 


^n,<7,(7' — ^n-\-‘2,(7,(7' ■) 


(S5) 
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since &aj^b,j Pn^N-n IV’)^ ® 7^ 0. Such a relation uniquely defines any ground state for any fixed local parity 

(even-even, odd-odd, even-odd, odd-even), and thus, for each fixed particle number N, there is a double degeneracy 
related to distinct wire parities. For example, a general ground state for 2N particles is given by 

\M (X 5] [woPLPkN-n) \i^): ® \K + WlPL+lP2N-(2n+l) Wo ® IV')^ , (S6) 

n 

and is parametrized by the complex coefficients wq and rci. This writing states explicitly that \4>2n) of the states 
\^L{N))ee \'^l{N))oo prosented in the main text. 

An alternative viewpoint on the states I^at) stems from considering two decoupled Kitaev chains at half filling 
/r = 0 and A = J. After defining the state |G) = |'!/')g ® where |'!/')g is the even-parity ground state of the wire 
a, the four ground states of the model read: 

{\G) \G) \G) \G)}, (S7) 

and are related to the edge Majorana fermions = ^ 2 l + *7?> where {%“,7^} = 2SijSafi, l 2 j-i = ~ “])> and 

72i = Ofj -|- dj, for a = a, b. It is possible to explicitly verify that the zero-energy subspaces of the number-conserving 
Hamiltonian H\=i are given by 

span { Pw IG), IG)} for N even, 

& span {Pat^^ IG) , \ G) } for N odd. 

Periodic boundary conditions 

The situation changes dramatically when one considers periodic boundary conditions. Let us consider for simplicity 
the case of an even number of sites. Then, for even values of N there is one zero-energy state, which is a ground 
state of the Hamiltonian. For odd values of N, no zero-energy state appears even in the thermodynamic limit; this is 
consistent with the fact discussed in the main text that the single-particle correlators decay exponentially, signaling 
the presence of a thermodynamic gap for the addition of a single particle. 

It is a remarkable and non-generic property of this model that the ground state for N even can be formally written 
in the same way as |V’i(-^))oo ^7- (4) main text. The demonstration is based on the fact that (see 

Eq. (1) in the main text) is the unique ground state of a Kitaev chain with periodic boundary conditions (L even) 
and follows the lines of the demonstration for hard-wall boundaries. Again, it is interesting to observe that: 

\MN))oo^PnWo®Wo (S8) 

even with periodic boundary conditions. The absence of an even-parity ground state for the periodic Kitaev wire 
provides intuition of why the odd N sector of the Hilbert space is gapped. 

ENTANGLEMENT SPECTRUM 

In this section we provide the detailed derivation for the entanglement spectrum of |V’L(-^))ee presented in the 
main text. We consider the reduced state of I sites on each wire pi = Tr(A_q [ |V’L(-^))ee (^i(-^)lee ] following 

expression identity operators on the first I sites are omitted): 

Pi = {jm\a ® (9m'It) IV'i(-^))ee (V’i (-^) I ee lim)a ® Wm') f, i (S9) 

where {jm} and {^} represent all possible fermionic configurations in L — i sites. Let us now explicitly compute 
{jm\a ® (9m'It, IpLiWeei ^^r even values of Ni = N -m-m': 

(J2m(2m+l)|a 0 (9tV-tV, -2m(2m+l) If, \MN))ee = jt/^t (iVt)),,(„„) i (SIO) 




Figure SI. Entanglement entropy Si for a block with I sites, {upper panel) considering a fixed filling n = ^ = ^, {lower panel) 
or a fixed particle number sector N = 10. 


which remarkably does not depend on the specific j or q. A similar relation exists for Ni odd, 


AC 


{32m(2m+l)\a ® {^N—Ni—2m(2m+l)\b (A^))ee l/^ 


ee,l,Ni 


ee,L,N 


\MNi)), 


3(oe) 


(Sll) 


Summing up such terms, we obtain the reduced state in diagonal form and its eigenvalues, presented in the main text. 

The demonstration for the double degeneracy in the entanglement spectrum in the limit of large L and L — I (i.e., 
large lattices and bipartitions not close to its edges), is related to the fact that, in this limit, Mee,i,Ni ~ J^oo,i,Ni- Even 
if we do not have an explicit analytical proof of the previous relation, numerical tests in several regimes corroborate 
this intuitive result. 

From the eigenvalues computed in this section, we can also compute the entanglement entropy of the block matrices, 
and see how it scales with the size of the block. We see in Fig. [^a behavior typical of a gapless Hamiltonian, which 
does not scale as an area law. 


EDGE MODES 

As discussed in the main text, in order to directly characterize the localization length of the edge modes we 
compute the following matrix element: Oj = ('!/'l(A^)Ioo IV'L(A^))ee consider here for simplicity the case of an 
even number of particles). This task reduces to a particular counting of the suitable fermionic configurations, starting 
from the explicit expressions in Eq. (4). 

If we act with the Vj = d^S_, operator on \ipLiN))^^, the only states 

which are not annnihilated bt Vj are those which have a particle at the j-th site of wire b and a hole at the j-th site 
of wire a. After the action of Vj, each of these configurations is changed by moving the fermion from the j-th site of 
wire b to that of wire a. Due to the anticommutation properties of fermionic operators, the state obtains the phase 

b ^ 

(—l)("«+"’i-), where ^ n“ is the number of particles located at the right of the j-th site of the wire a, and 

r=j + l 

f-1 

the number of particles located at the left of the j-th site of the wire b. These phases describe the parity 

r—1 

of the configuration on the segments [j -|- 1,L] for the wire a and [l,j — 1] for the wire b. Since the ground state is 
an equal superposition of all configurations distributing N particles between the two wires and fixing the parity of 


j 2 n) 0 |gAr- 2 n)a (see Eq. in the main text) 
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the number of particles for each wire, varies from a minimum value equal to max(0,-/V — j) to a maximum value 
min(A^, L — j)- Analogous relations exist for The matrix element Oj is thus related to the simple counting of 
particle configurations taking into account the phase (— 

Let us discuss this in more detail. The total number of fermionic configurations which originate the state 
is 



b—wire a—wire 


(S12) 


whereas those which have a particle at the j-th site of wire b and a hole at the j-th site of wire a are 


E 


L- 1\ / L- 1 
2n-l) [N-2n 


= E E 


j -l\ f L- j 
2n — 1 — 


E 


L- j 


3 - 1 


J \N — 2n — 


b—wire a—wire 


b—'L 


a—wire 


(S13) 


Taking into account the phase (—we have that 




\/-N'ee,L,N-N'oo,L,N u ^ 
'' n,n^,n^ 




L-j\ f j - 1 
N — 2n — 


b—wire 


a—wire 


(S14) 

Following the same steps, it is not hard to see that for a more general number-conserving two-body operator we 
have: 


{MN)loa]br\MN))ee = 


\/-l^ee,L,N-N'c 


E (- 1 )“ 




L,NJ^oo,L,N u „ 


r — 1 


L — r 
2n — 1 — 


L- j 


J - 1 

N — 2n — n' 


R. 


b—wire 


(S15) 

for (L -|- r) — j > 1. In particular, edge-edge correlations have a particularly simple expression in the thermodynamic 
limit. Let us demonstrate that starting from: 


{MN)\ooblbL\Mm)ee 

{MN)\ooblh\MN)),, 


1 


\/-^ee,L,N-^oo,L,N 

1 


E(-i) 


(JV-l) 


E 


L- 1 



y/j^ee,L,NJ^ooX,N „ — 1/^ \N — 2n)^ 

b—wire 


a—wire 


(516) 

(517) 


Using the Chu-Vandermonde identity, X]/c=o (T) (r-fe) “ which holds for non-negative integer m,n,r, we 

obtain that in the limit of large lattices the following is true: 


" '-V--V-^ 

b—wire a—wire 


(S18) 


Thus, edge-edge correlations reduce to: 

{MN)\oob\bL\MN)), 


Similar expressions exist for {'tf^L{^)\oo I^L(-^))ee- Note that, if N is odd^ a minus sign appears due to the overall 

phase in Eq. (|S16D: {ipLiN)\^^a\bL |^/’L(IV))ee 


v{\ - v) L_i.oo 

, if. 

4 -v{l-v). 




(S19) 

























10 


SINGLE-PARTICLE AND SUPERFLUID CORRELATIONS 


In a general way, any ground state observable can be computed as in the previous section through a simple counting 
of suitable configurations. In this section we evaluate the single particle correlations {Ojdj+r), as well as the superfluid 
correlations {d\d\_^-y^dj+idj). We skip unnecessary details and focus mainly on the presentation of the final results. 
We only consider the ground states for even values of N because the odd case is mathematically equivalent. 


Single particle correlations: 


(V'L(^)lee(oo)«i%+>- I V'i (^))ee(oo) = 


(S20) 




(oo),L,N 


E (-1)' 


(j.r) 


(j.r-) 


L 


r — 1 


L — r — 1 


'V2n(2n+1); \N - 2n{2n + 1) - 1 - 


b—wire 


-V— 

a—wire 




■E 


f ^ 


(oo).r,« ^ \2n{2n + 1)J \N - 2n(2n + 1) - 1 


L — r — 1 


if r > 1 ; 


if r < 1 . 


(S21) 


b—wire 


a—wire 


Here, = Ei=i number of particles between the sites j and j + r of the wire a, which varies from a 

minimum of zero (where all the particles are in the wire 6 ), to a maximum value equal to min(fV — l,r — 1) (where 
either all the remaining N — 1 particles lie between these sites or the wire segment is completely filled). 


P-wave superfluid correlations: 


Pairing correlations characterize the model defined by Eq. (3) and are of p-wave nature. Their behaviour is captured 
by the pairing operator dj^idj. At the exactly solvable point, A = 1, these correlations can be computed analytically. 
The evaluation proceeds as in the previous cases, relying on the counting of the suitable fermionic configurations. The 
result for the p-wave correlation reads: 


(V’i.(^)L(oo) 44+iaj+iaj |V’i.(A^))ee(oo) _/v;e(oo)(A, A^) ^ - 2n - 2) 


i + l<j. (S22) 


b- 


wire a—wire 


Using the Chu-Vandermonde identity, we obtain in the limit of large lattices, 

' L\ / L - 4 \ 1 /2L - 4 


E 


2nJ \N-2n-2j 2\N-2 


b—wire a—wire 


with A/*ee(oo) (-^ 7 -^) ~ Mat)’ large distances the p-wave pairing correlations saturate to 

IV'L(A^))ee(oo) ~ -^ (1 - 0 • 


242.4 + 0 ( 2 , 3 ) 


(S23) 


(S24) 


This non-zero value for large \i—j \ distances hints to the fact that there are no sound-like modes in the exactly-solvable 
model, since they would lead to the algebraic decay of the p-wave correlations. 
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BRAIDING 


In this section we demonstrate the properties of the braiding operator BaR,bR{j) = (^ + ZaR^tRj )/presented in 
the main text. To that end, it is useful to first study the action of ZaR^bR,j, starting from the simpler case j = 1: 


ZaR,bR,l = ; XaR^bR.l = o'r^L - -^1 = (^^ 5 ) 

J~l ZL/ 

The action of a\^bR on \ipLiX))^^ produces an unnormalized state proportional to the equal-weighted superposition 
of all fermionic configurations with (i) N fermions, (ii) odd wire parities, (iii) a hole in the L-th site of the wire a, 
and (iv) a particle in the L-th site of the wire b. Due to the anticommuting properties of fermionic operators, this 
state gets a global phase (—l)(^“i)+(^“) = (—(recall that (—1)'^“ and (—1)'^'’ are well-defined although Na 
and Nb are not fixed). With similar reasoning one can also characterize b\aL (here the hole (fermion) is in 

the L-th site of the wire a (6) and the global phase is (—l)('^““i)+(^“i) = (—l)^*). The state ZaR^bR,i |V'L(-^))ee 
the normalized state described by the equal weighted superposition of all fermionic configurations with (i) N fermions 
and (ii) odd wire parities and without (iii) the simultaneous presence of two holes or two particles at the L-th site of 
both wires a and b. The normalization constant Fi corresponds to the “fidelity” of the state with the ground state 

with such local parities, (iAl(A^)Ioo ZaR,bR,i \'^l{N))^^ = y/l - + (1 - 

The operator ZaR,bR,j acts on the last j sites of the wires. Let us consider for example j = 2: 


-'aR,bR,2 — 


^aR,bR.l “ 1 “ ^aR,bR.l^aR.bR,2 ^ 
^2 ’ 


YaR,bR,lXaR,bR,2 = {nln^R + (1 “ 0(1 “ ^^L)) (oL-I^L-I “ b^R_^aL-l^ ; 


F 2 = v'l- [j^O(l-J^)4]. (S26) 


Note that the term YaR.bR^iXaR^bR ,2 \'4’l{X ))acts only on the fermionic configurations which were missing in 
XaR.bR,! \'^l{N)) However, it is clear that ZaR^bR ,2 |'*/’L(-^))ee does not contain any configuration with four holes 
or four particles in the L-th and (L — l)-th sites of both wires. Whereas this still makes ZaR,bR ,2 \i^L{N)),, dif¬ 
ferent from |V'i(N^))oo, it is a considerable improvement with respect to the previous case. In general, the state 
ZaR,bR,j |'0L(-^))ee equal weighted superposition of all fermionic configurations with (i) N fermions and 

(ii) odd wire parities, and without (iii) the simultaneous presence of 2j holes or 2j particles in the last j sites 
of both wires. L) corresponds to the fidelity of the state to the ground state with such local parity pattern, 

j |'0L(lV))gg = -I- (1 — iz)2t]. For a large enough j < L/2, the difference between 

ZaR,bR.j |'0L(lV))gg and |V'i(-^))oo becomes exponentially small, so that the above operator in the zero-energy subspace 
reads 


ZaR,bR,, ~ E i\MN))ee {MN)lo ” {MN)U + 

N even 

+ i\MN))oe {MN)\eo - \MN))eo {MN)U ■ (^27) 

N odd 


For j and L are such that the error is negligible, we define the number conserving operator RaR,bR = BaR^bR{j)- 
We observe that RaR.bR acts on the zero-energy states in the same way done in the conventional number non¬ 
conserving scenario of two neighboring Kitaev wires by the braiding operator RaR bR — l^RlhR — {I->r^aRlbR)/'/^, 
where are the zero-energy Majorana operators exponentially localized at the A = R, L edge of the wire m = a,b. 
In order to verify this explicitly, we first recall that the number non-conserving edges Majoranas are related to a 
non-local fermion as fm = jmL — ilmR which is the Bogoliubov zero-energy mode. The two degenerate ground states 
of the wire m = a, b, \tp).auT ^be parity of the number of fermions, even, e or odd, o, and labels the two ground 
states), correspond to the presence or absence of the non-local fermion /„: /m \'4’)rae ~ b, and IV')me = IV’)mo' 

Using the inverse relations y-mi oc ^/m -I- /m) and ymi? oc i (^fm — /m)) b is now easy to see that: 

%RlbR |V')a<T IV’)br = P< 7 {%R \^)ao)i%R \^) br) = (-*P^) (“fPr ) |V’)aff \^) bf = -Pr IV')aff \^) br (^ 28 ) 


where Po- = 1 for a = e and pa- = —1 for a = o\ a, t are the flipped tr, r. As claimed, jaR%R acts on the ground space 
in a way which is completely analogous to that of in Eq. (S271 for the number-conserving model. The equivalence 
of RaR.bR and R'aR^hR foUows directly. 
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The unitarity of the braiding operator RaR.bR restricted to the ground subspace can be proved explicitly. Let us 
first notice that: 


{i’L{N)lrKR,bRRaR,bRHL{N)),.^, (IV) (l + . 


^aR,bR,j 

t 

aR,bR,j ^ 


aR,bR^j 


“'aR,bR,j 


jZaR,bR,j] \lpL{^))cr>r^ 


^aR,bR,jj — ^(T,a''^T,r 


(S29) 


where in the second line we used the fact that ZaR^bR,j is anti-Hermitian. Thus, in the ground state subspace we 
have, 

^9 ^aR,bR RaR,bR PI = Pg (S30) 

With similar procedures we can define other braiding operators with completely analogous properties. For example, 
let us consider the transformation z o (L +1 — z) which maps an operator at site z to site (T +1 — z), thus mapping the 
right edge to the left one and viceversa. When we apply it to ZaR^bR,j, we can define ZaL,bL,j, which is exponentially 
localized at the left edge of the wire. For j and L large enough, its explicit expression is: 

ZaL,bL,j^ E i-\MN))ee{MN)\oo+\MN))ao{MN)\ee) 

N even 

+ Y. i\MN))oe {MN)\aa ” \MN)) {MN)\oe) ■ (S31) 

N odd 


The braiding operator RaL,bL can be defined, and an explicit calculation shows that it is unitary and that it resembles 
the operator R'^j^ for the number non-conserving case. 

As a second example, let us consider the transformation —>■ —iaj for j < L/2 which maps the right edge of 

the wire b to the left edge of the wire a, leaving the other two edges unchanged. When we apply it to ZaR^bRjt we 
can define the operator ZaR^aLj, whose explicit expression for j and L large enough is: 

ZaR,aL,g ^ Y * (\MN))ae {MN)\ee " \MN))oo {MN)U 
N even 

+ ^ Z (|V'L(iV))eo (V’L(^)leo ” I^L(^))oe {^^ 00 ) ■ (S32) 

N odd 


Again, everything follows as before. In general, with this method one can define sixteen braiding operators RmA,m'A, 
with m, m' = a, b labeling the wires and A, A' = L,R labeling the left and right edge. 

Let us conclude with an explicit verification of the non-abelian character of these operators; to this aim, we initialize 
the system in the state |V'L(-^))ee- then perform two braiding operations on the edges in different sequences: 

(^RaR,aLRaR,bR — RaR,bRRaR,aL^ |V’L(-^))ee = 2 ZaR,bR,j — ZaR^bR,]ZaR,aL,i^ |V'L(-^))ee = * IV'l(-^))oo ■ 

(S33) 


The result shows explicitly the non-commutativity of the braiding operations. 



